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Abstract 

In a recent paper Stevenson claimed that analysis of the data on the wave function 
renormalization constant near the critical point of the Ad Ising model is not consistent 
with analytical expectations. Here we present data with improved statistics and 
show that the results are indeed consistent with conventional wisdom once one takes 
into account the uncertainty of lattice artifacts in the analytical computations. 
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1 Introduction 



One of the apparently simplest quantum field theories in four dimensions is the <p 4 
theory with n components. Conventional wisdom (CW) holds that the theory is 
trivial for all n > 1 . Unfortunately there is presently no rigorous proof. One 
relies heavily on the validity of renormalization group (RG) equations for physical 
quantities [T] together with boundary conditions at finite cutoff provided by non- 
perturbative methods |2|3j . 

Apart from its purely theoretical interest it is of phenomenological relevance 
since for the case n = 4 it constitutes the pure Higgs sector of the Minimal Standard 
Model (MSM). The fact that cj>\ is (probably) trivial does not invalidate (renormal- 
ized) perturbative computations for amplitudes at energies well below the physical 
cutoff where the MSM may be a good effective theory. 

In the past triviality has been invoked to propose upper bounds on the mass 
of the Higgs boson (see e.g. refs. (IHE1)- It must be stressed that these bounds 
are non-universal, they depend on the particular regularization. But for a given 
regularization it is conventionally accepted that such a bound can be given a pre- 
cise meaning. In recent papers Cea, Cosmai, and Consoli (CCC) |9|10| claim that 
triviality itself cannot be used to place upper bounds on the Higgs mass even for a 
given regularization. They assert that standard predictions of the RG analysis for 
the behavior of some quantities near the critical line are not valid. If true this would 
indeed be rather important because it would reveal a serious flaw in our conventional 
theoretical understanding of the pure </>|! 

In ref. Duncan, Willey and the present authors explained why critiques of 
the standard picture raised in ref. |5] were not relevant. Nevertheless one must admit 
that the unconventional picture of CCC cannot be ruled out by present numerical 
simulations. Recently two papers appeared, the first by CCC and the second 
by Stevenson [T5|, again claiming finer but significant discrepancies between quan- 
titative (standard) analytic predictions and numerical data in the 4d Ising model 2 . 

It is the purpose of this paper to reply to these challenges and to demonstrate 
that they are too weak to seriously cast doubt on CW. The main objection by CCC 
and one by Stevenson are rather easy to dismiss. A second objection by Stevenson 
is more difficult. It concerns a certain difference A between the wave function 
renormalization constants below and above the critical point 

A = Z^{k = 0.074) -Z r {k = 0.0751) , (1.1) 

1 in the sense that if the model is defined non-perturbatively with an ultraviolet cutoff A, say via 
lattice regularization, connected r— point functions with r > 2 vanish in the limit A — » oo 

2 The Ising model is obtained as the limit of the 4> 4 model when the bare coupling goes to infinity. 
As such one generally considers this the "worst case" i.e. if CW holds for the Ising model it is even 
more plausible for finite bare coupling. 
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which we have called "Stevenson's step" in the title. The k values chosen here 
are about the closest to the critical point where one can presently obtain good 
statistics with some (but not unreasonable) computational effort 3 . Stevenson claims 
that "theoretical predictions for A cannot be pushed above 0.05 well short of the 
"experimental" value Amc = 0.071(6)". It is of course debatable whether such a 
small discrepancy indicates a potential problem, however we decided it merited more 
careful investigation. 

It is however clear that one is here addressing few percent effects, and to clarify 
the situation we need both data and analyses which are precise to this level. The 
main analytic sources of error concern the treatment of higher order cutoff effects in 
the framework of renormalized PT. The main sources of error in the numerical side 
are the determinations of the zero momentum mass m R ; apart from the statistical 
errors one has the systematic errors in the procedures to extract m R from the (finite 
volume) data particularly in the symmetry broken phase. 

In the next section we present a summary of the available raw data in both the 
symmetric and broken phases. We have performed simulations in both phases and 
in particular increased the statistics at previously measured k values in the broken 
phase by a factor of ~ 10. 

We then discuss various determinations of Z R from the data and compare with 
theoretical expectations. Next we show that these data are not in contradiction 
with conventional wisdom. The reason for this conclusion differing from that of 
Stevenson has two main origins. Firstly unfortunately our central value of Z R at 
k = 0.0751 in is about one standard deviation lower than that obtained from 
the present run. In this connection we remark that in those runs we were not aiming 
at high precision but only sufficient to reach our goal to present evidence that Z R 
is not increasing logarithmically as one approaches n c . Secondly we point out that 
there is a quantitative uncertainty on the 0(a 2 ) lattice artifacts which are of course 
increasingly relevant as one goes away from the critical point. This is of course not 
at all new, however sometimes forgotten and perhaps underestimated in standard 
RG analyses. 

2 Ising MC simulation and results 

We work on hypercubic lattices of volume L 4 with periodic boundary conditions in 
each direction and with standard action. In this paper we adopt the notations and 
definitions in ref. and will generally not repeat them here. 

If one just wanted to obtain Stevenson's step, measurements at only two k 
values are required. However precise simulations at these points are CPU expensive 

3 the correlation lengths are around 6 in lattice units. 
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and hence it is useful to compute also at points with smaller correlation length 
to observe the approach to Stevenson's chosen points which hopefully also reflects 
the approach to the critical point. As mentioned above we have simulated in both 
phases. In the symmetric phase this served as a check with respect to previous 
simulation results by other authors. In Table ^ we collect the (to our knowledge) 
best data for observables in the symmetric phase, which are obtained without a 
fitting procedure; by "best" we mean data with the largest statistics and for lattices 
with large physical volumes mL > 6. (Our tables can be found in Appendix B.) All 
the values in this table are from the present simulation except those for k = 0.07102 
which come from ref . |14j . The previous results of Montvay, Miinster and Wolff ^H] 
for their lattices A,C (k = 0.071, L = 12 and k = 0.0732, L = 20) are in complete 
agreement with ours, but have larger errors. 4 Note that we have also measured at 
a new point k = 0.07436 slightly closer to n c than previous ones. 
Also included in Tabled is the quantity 



mS aive 



ko 2 G{k ) , /2tt 



\ G(0) - G(k ) ' V L 



0,0,0 , (2.1) 



. naive 



where G is the Fourier transform of the (connected) two-point function, 
tends to the desired zero momentum mass in the infinite volume limit 5 . 

In Table |21 we present a similar data collection for the broken phase. This is an 
update of Table 3 in here we have increased our statistics by a factor of ~ 10. 
In these runs we also measured the connected 3-point function X3 (except for one 
lattice). For a large subset of the runs in the broken phase we also measured the 
correlation matrix of the time slice field cp(t) with the composite field 



<p 2 (t) oc^(t) 2 , = x). (2.2) 



3 Determination of m R and derived quantities 

The usual expressions for the wave function renormalization constant and renor- 
malized couplings involve the infinite volume zero momentum mass. In particular 
a precise determination of Z n (as required to discuss Stevenson's step) requires an 
equally reliable determination of m R . Firstly in considering the r-point functions one 
can verify using the formulae given in ref. ^Z] that finite volume effects coming from 
tunneling are negligible for our lattices. We then adopted two fitting procedures. 

4 Some of the entries for their lattice B (k = 0.0724, L = 16) are many standard deviations away 
from ours and almost certainly wrong (as also suspected by Stevenson \12]). 
and is volume independent for the free lattice theory with standard action 
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In the first we made fits of G(k)^ 1 for k = (2tts/L, 0, 0, 0)) for small values of s 
(s < 3 in the symmetric phase and s < 4 in the broken phase) to polynomials in k 2 
and just recorded good fits with x 2 /dof < 1. 

Due to the discreteness of available values of k the determination of the slope 
at k = is prone to discretization error. A removal of this deficit is attempted in a 
second procedure where we performed fits of the connected (time slice) correlation 
function S(t) to obtain the physical mass m. Relying on the expectation that finite 
volume effects on m are extremely small in this model when mL > 7, we obtain an 
estimate for the infinite volume m R by computing the second moment with the data 
in a large portion of the lattice volume 6 and then computing the contribution from 
the rest of the (infinite) volume using the measurement of m. 

Our fits of S(t) were in terms of one and two mass cosh functions; for the 
1— mass case we fitted only distances t > 1/m and for the 2— mass case t > 2/3m. 
The 2— mass fits were constrained to have the second mass fixed to 3m for the 
symmetric phase and 2m for the broken phase. The 2— mass fits give (as expected) 
a slightly lower central value of m but a slightly larger error than the 1— mass fit. 
This is again reflected in the resulting values for m R and Z R . Figs. ^ El are typical 
examples in the symmetric and broken phases respectively, which illustrate the good 
quality of the fits and the relatively very small contribution of the higher particle 
states. 

For the subset of data in the broken phase where we had the full correlation 
matrix mentioned above we found that the two operators cp, ip% were nearly parallel, 
both coupling very weakly to the 2-particle state, and thus this did not help to 
significantly reduce the errors on m. 

In Tables El and 0] we give results for m and quantities derived from various 
estimates of m R in the symmetric and broken phase respectively. In practically all 
cases all methods to determine m R gave compatible results 7 . 

In Fig. |2l we plot renormalized couplings (a R = <7 R /(16-7r 2 )) in the symmetric 
and broken phases using m R = m H aive versus 

2sign(r) k 

«oW = Q1 1 1 , T = 1 , I 3 - 1 ) 

3m|r| k c 

with the presently best estimate of k c = 0.074848 from ref. ^H]. In the symmetric 
phase the coupling is defined through the connected 4-point function. In the broken 
phase we have included the renormalized coupling g R defined through the vacuum 
expectation value and another coupling defined through the connected 3-point 



6 e.g. t < 3/m where the errors are reasonably small 

7 One exception is the lattice k = 0.0752 where a best fit gave a slightly higher value of m R . 
The resulting value of Z R is bigger than that at k = 0.0751 which we consider as a signal of the 
potential instability of such momentum space fits. 
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Figure 1: Ratio of the measured S(t) to the 1— and 2— mass contributions (filled circles and 
squares respectively) to the 2— mass fit for the k — 0.0744, L = 52 lattice. 




Figure 2: As for Fig.[T]but for the n = 0.0751, L = 48 lattice. 
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Figure 3: Renormalized couplings versus ao; ao < corresponds to the symmetric phase and 
ao > to the broken phase. Full circles and squares represent qr and Qr3, respectively. For 
ao < we have plotted — qr. 

function. According to standard RG analysis these should behave as 167T 2 |cto | as 
r — > 0. Note that the relations of the measured values for g R and g R s are quite 
consistent with those obtained in renormalized perturbation theory: 

Sr3 = 9r (1 + 3a R + 3.75a R + ...). (3.2) 
3.1 Stevenson's step 

Stevenson presents plots with respect to k whereas we prefer to present plots wrt 
ao since dependence of quantities of interest on this variable are expected to be 
smoother. In Fig. 21 we plot Z R obtained using m R values obtained from the 1-mass 
fit method. Using m R = m R aive or m R corresponding to the 2-mass fit method 
would result in very similar plots, with somewhat larger errors for the latter case. 
(The same applies also to Fig. [SJ) Here one certainly does not see any signal of a 
discontinuity at the critical point. As mentioned above the value of Z r (k = 0.0751) 
cited in is over one standard deviation away from our present value 8 . For 
Stevenson's step we would now quote a measured value of 

A MC = 0.057(6) . (3.3) 

8 Also note that our present value of Z R at k = 0.0754 has moved down wrt to that quoted in 
and is now consistent with some as yet unpublished data by Cea and Cosmai mentioned in 



6 



"i 1 r 



"i 1 r 



0.95 ■ 



0.9 ■ 



J i L 



i 



J i L 



-0.2 -0.1 0.1 0.2 



Figure 4: Measured values of Zn using raR obtained by the 1-mass fit method. 



Given the more precise value of Z R at k, = 0.0751 Stevenson may not have written 
his paper. On the other hand this new value for A is still bigger than what Steven- 
son calls analytically feasible. We however disagree with this opinion. The problem 
is that our quantitative control of 0(a 2 ) effects in the RG equations is not sufficient. 
Although these are rather small at the values of k used in the definition of A, they 
are still not negligible when discussing possible discrepancies between theory and 
"experiment". In refs. |2J3| some such effects were taken into account by includ- 
ing the m R dependence of perturbative RG coefficients appearing at low orders of 
perturbation theory. This procedure was also adopted by Stevenson \12\ . It is how- 
ever just a pragmatic procedure (i.e. practically the best one could quantitatively 
do at the time), but it is not a quantitatively systematic prescription. Firstly it is 
not consistent to include 0(a 2 ) effects while ignoring higher perturbative effects 9 . 
Secondly even if one disregards this, the leading 0(m|) effects at r-loop order are 
of generically of the form g R +r m R \ hi?7iR.| r and hence all quantitatively of the same 
order since the RG equations predict g R ~ l/lnm R . 

To obtain the leading cutoff effects one must follow the method of Symanzik ^3] . 
The analysis shows that the leading artifacts for the correlation functions (for the 
case n = 1) are of the form 

r ~ T (0) +0 (m R |lnm R | 1/3 ) . (3.4) 

9 analogous procedures are often (similarly questionably) adopted in phenomenology when taking 
higher twist effects into account. 
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Here is the formal perturbative sum; in the process of defining this non- 
perturbatively it could be that renormalon-type effects would lead to cutoff effects 
of the same form as the leading operator insertion. But even though the form of the 
leading cutoff corrections may be known the amplitude is undetermined. 

Another simple exercise to appreciate this point is to compute in the leading 
order of the 1/n expansion. In that limit the cutoff effects are of the order m R 
however these are not given by taking the limit n — > oo of the first perturbative 
contributions. Some illustrations are given in Appendix A. 

Having realized that we have unfortunately insufficient quantitative knowledge 
of m^-effects it is still legitimate to ask if the numerically found value of A looks 
inconsistent with the conventional theoretical expectation that Z R approaches the 
same value C 2 coming from both sides of the critical point: 

CW : lim Z~ R = C 2 = lim_ Z~ R , C 2 = 2k c C 2 . (3.5) 

To illustrate the situation we consider two expressions 6 { 2 I] and (5f 7) which we 
would eventually expect to approach C 2 faster than Z R . The first, which arises 
in the renormalization scheme of is defined by dividing Z R by its perturbative 
expansion truncated at 2-loops: 



Z R ( l + -^a R + 0.100896a| ) , for k < n c , 



1 

18' 
7 

36' 



-1 



Z R [1- —a R - 0.538874a; ) , for k > k c . 



-1 



(3.6) 



The second, which is a natural choice in a field theoretical context, is merely defined 
from the first by 

di H) ee AgW . (3.7) 



c 



The difference between the two functions is just an order m R cutoff effect |2I3| : 

k c - k = C 3 m R g R 1/3 {1 + 0(g R )} , for k < n c , (3.8) 

k-k c = -C' 3 m^g R 1/3 {1 + 0(g n )} , for k > k c , (3.9) 

which incidentally has the same cutoff effects as in Ij3.4|) . In Fig. [5] we plot them 
together to give an idea on the importance of the 0(m R ) effects. Both are certainly 
not inconsistent with the expectation that the limits are the same on both sides 
(note the difference in scales on the vertical axis) . Indeed if one allowed to naively 
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Figure 5: Top, c!jp (renormalization scheme of refs. [2,3]). Bottom, C^ 1 ^ (field theoretical 
renormalization scheme). Note the different vertical scales. 
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extrapolate the curves by eye one would probably infer a different sign for the dif- 
ference of the limiting values from the two plots. 

If one had much more precise data one could attempt constrained fits to the 
CW but this is not warranted with the present data. Assuming CW we would now 
quote C 2 = 0.95(1) which corresponds to lnC2 = 1.85(1) to be compared with the 
result lnC 2 = 1.87(1) in 0. 

4 Reply to some other critiques 

In ref. ^Hj the authors reconsider the quantity v 2 x i n the broken phase which ac- 
cording to CW behaves as 

v 2 X 

where £ = | ln(/c — k c )\ with 

«i 
a 1 

A fit of the expression flUTJ) to the data gave [H] oi = 1.267(14) and a 2 = -2.89(8) 
whereas the theoretical prediction based on the results quoted in jHj for the values 
of the non-perturbative constants C\, C 2 and C3 is a\ = 1.20(3) and a 2 = —1.6(5). 
The authors of ^0] claim that such a comparison "shows that the quality of the 
2-loop fit is poor". How they can reach such a conclusion is surprising to us. Firstly 
the values a\ agree within one standard deviation. Secondly the value of a 2 from 
the fits can only be regarded as effective. Higher order terms e.g. of the form 1 jl 
are completely neglected in the fit. 

We have made new fits including all data from Table EJ omitting only the 
K > 0.076 lattices and the one corresponding to k = 0.07504 (because of its too 
small physical volume). Our first fit, which is shown in Fig. El has an acceptable 
X 2 /dof = 1.2 and yields 

a x = 1.224(6) and a 2 = -2.68(3). (4.4) 

As discussed at the end of the previous section, our Z n data prefer a slightly smaller 
C 2 value corresponding to lnC 2 = 1.85(1). This changes the theoretical prediction 
to ai = 1.15(3). We have made a second two-parameter fit with a third term of the 
form as/l added but where the first coefficient was kept fixed at a\ = 1.15. The 



a (£-— hit) +„■> + . 



(4.1) 



9C 2 2 

327T 2 ' 

InCg + 21nC{ - 1.6317. 



(4.2) 
(4.3) 
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Figure 6: 2-parameter fit of the expression 1)4. Q to the data from Table El 

result of this fit 10 is 

oi = 1.15, a 2 = -1.80(3), a 3 = -3.5(2), xV^of = 1.0. (4.5) 

The fact that it is impossible to distinguish between (|4.4f) and (|4.5|) illustrates the 
point we made in the previous paragraph. Taking into account all the uncertainties 
we find the agreement between theory and MC 'experiment' satisfactory. 

We completely disagree with the implications of Section 5 in ref. ^2]- Here 
Stevenson proposes that a proper extraction of m R is obtained by globally fitting 
the data for G{k)~ l in the whole available momentum range to the form of the 
free lattice propagator with standard (nearest neighbor) action. But this is an a 
priori incorrect procedure since m R must be extracted from data including only 
low momenta. Accepting this fact it is then a rather surprising "experimental" 
finding that the numerically measured propagator is so extremely close to the naive 
propagator; only a detailed look reveals that there is some significant deviation at 
larger momenta. It is a bit easier to see the deviation in coordinate space e.g. in 
Fig. 0we plot the effective masses at ft = 0.0751. The effective mass m e g (t + 0.5) is 
defined for the 1-mass case by mi in the ansatz B [exp(— m\t') + exp(— m\{L — t'))] 
with the parameters defined from S(t') at t' = t,t + 1. This would be constant for 
a free standard lattice propagator, but the data shows clear i-dependence. For the 
constrained 2-mass fit Y2i=i [exp(— mrf') + exp(— m,(L — t'))] with m<i = 2mi, 
m cs(t) = fn\ from the correlation function at if = t — 1, t, t + 1. 
10 which has incidentally a value of a,2 quite close to the prediction of ref. 
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Figure 7: Effective masses for the k = 0.0751, L = 48 lattice. The circles are from the 1-mass 
definition, and the squares from the 2-mass fit definition as described in the text. 



5 Conclusions 

An experimental observation contradicting a prediction of an until then accepted 
theory is always an exciting event. It invalidates the theory as it stands and in- 
evitably leads to progress in our understanding. Similarly finding mismatches be- 
tween theoretical predictions and numerical simulations in the cftf theory as claimed 
in refs. |1UI12| would be a serious blow if they withheld scrutiny. We hope to have 
convinced the reader in this paper that conventional wisdom concerning this struc- 
turally simple theory is still alive. Although present numerical simulations support 
CW, the scenario can unfortunately only be "nailed down" by analytic proofs. 
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Appendix A: leading order 1/n expansion 

For the n-component model we take over the notations of ref. [7]. In the 1/n expan- 
sion one takes n — > oo with 

A = nA (A.l) 

held fixed. 

In the symmetric phase in leading order we have 

Zn = 1 , (A.2) 

and given A, k < k c = ^Ji(O) the renormalized zero momentum mass m H is deter- 
mined by 

« 2 ( 2 - 1 



y (m R + 8-- ) +2k- Ji(m H ) = 0, (A.3) 

where (see ref. j2]) 



d 4 fc / , o 2 



U (2*0 

The renormalized coupling (fo = ng R is given by 



6 

to = 73 — • (A.5) 

y + J 2 (m R ) 



For A = 00 there is a simplification and 



dfe _ 24m R J 3 (m R ) 

mR a^ " j 2 (m R )2 (A - 6) 

= y^ 2 + m R 5i(m R ), (A.7) 



with si = l/(16vr 2 ) and 

24 ( j 3 (m R ) - 5^) 
J 2 (m R 

With the expansions of J q (m R ) in [2] for small m R the non-scaling piece behaves as 



Bi(m 8 ) = V j ( ^ A2 • (A. 8) 



#i(m R )~ r ^-a T + ... (A.9) 

s\ m(m R ) 
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In renormalized PT one obtains (in the limit n — > oo): 



dm R 



Ami 



+ m\ 



-9n + 9r 



2 f si 2m R 



Ji(m E 



+ m 



2\2 



So including the leading non-scaling effects we write 

,2 



?R Si— 2 



+ 0(;^ 3 m R ) . 

(A.10) 

(A.ll) 



where f s denotes the coefficient from the s loops. For the first coefficients we have 

3 



/« 



o 



3 1 



fx 



_3_ 

si 



(A.12) 



to be compared to the non-perturbative result in Eq. (|A.9|> . 

Although Z R is so simple, the 7-function defined in [2j is not zero e.g. in the 
limit A — > 00: 



u = — m R — In Z R = m R 



"5m 



R 



(A = 00) . 



For small m n we have 



■^m|lnK) + ... 



(A.13) 



(A.14) 



whereas in renormalized perturbation theory one obtains in the leading order 1/n 
expansion 



mi 



8 + ml 



</l(m R ) 

8 + m R 



mi 



1 + 1 + 0(fe 2 



(A.15) 
(A.16) 



which has no m R ln(m R ) behavior at tree and 1-loop level whereas the full non- 
perturbative function does. 

Similar features are found in the symmetry broken phase. 
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Appendix B: Tables 



K 


L 


X 


~X4 


naive 
R 


0.071 


12 


27.9884(29) 


5.602(26) • 10 5 


0.49528(4) 


0.071 


18 


28.0184(12) 


5.606(26) • 10 5 


0.49500(2) 


0.07102 


12 


28.20(3) 


5.9(6) • 10 5 


0.4924(5)* 


0.0724 


24 


45.697(14) 


3.68(19) • 10 6 


0.38330(14) 


0.0732 


20 


69.991(21) 


1.749(44) . io 7 


0.30760(8) 


0.0732 


30 


70.009(9) 


1.768(50) • 10 7 


0.30764(4) 


0.074 


40 


142.676(62) 


2.74(23) • 10 8 


0.21392(9) 


0.07436 


52 


257.27(19) 


2.14(45) • 10 9 


0.15871(12) 



Table 1: Measured values of X>X4, and m^ alve in the symmetric phase of the Ising model. The 
star indicates that for this entry ma™ was not measured directly, but could be estimated from the 
published data. 
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K 


L 


X 


V 


naive 
R 


~X3 


ref. 




10 


5.130(2) 


0.571267(8) 


1.0041(5)+ 




HU 


0.077 


16 


18.196(9) 


0.389500(9) 


0.5595(5) 


3.734(22) 


10 3 




0.077 


32 


18.21(4) 


0.38951(1) 








m 


0.076 


20 


37.80(3) 


0.301544(13) 


0.3942(5) 


1.942(14) 


10 4 




0.076 


32 


37.70(31) 


0.3015(1) 


0.428(5)* 






M 


0.0759 


32 


41.71(13) 


0.29030(2) 










0.0759 


48 


41.95(93) 


0.29028(5) 








M 


0.075628 


48 


58.70(42) 


0.25580(2) 






M 


0.0754 


32 


87.08(13) 


0.220488(14) 


0.2622(6) 






0.075313 


48 


104.2(1.3) 


0.20477(4) 






M 


0.075231 


60 


130.8(1.4) 


0.18812(3) 






M 


0.0752 


36 


142.1(8) 


0.18138(5) 


0.2054(16)+ 






HU 


0.0752 


40 


143.03(28) 


0.181291(14) 


0.2054(6) 


4.374(80) 


10 5 




0.0752 


48 


142.6(9) 


0.18132(4) 


0.2055(18)+ 






mi 


0.07512 


32 


193.1(1.7) 


0.1617(1) 


0.206(4)* 




m 


0.0751 


48 


206.32(40) 


0.156532(12) 


0.1715(5) 


1.043(19) 


10 6 




0.0751 


52 


201.2(6.2) 


0.15654(7) 








M 


0.0751 


60 


202.4(8.6) 


0.15648(2) 








M 


0.07504 


32 


293.4(2.9) 


0.13822(12) 


0.172(3)* 




m 


0.074968 


68 


460.2(4.9) 


0.11261(5) 








0.0749 


68 


1125(36) 


0.07736(12) 








m 


0.0749 


72 


1141(39) 


0.07752(21) 











Table 2: Measured values of x> u > m R alve an< i X3 from various Ising simulations; data from this 
investigation have no entry in the last column, kq — 0.080795. In the mR column + indicates that 
instead of the naive, the 1-mass fit method has been used and * indicates the mi a tt mass of ref. [16]. 
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K 


L 


iter 


m 


m R 




9r 


0.071 


18 


34M 


0.49007(2) 


0.49499(2) 
0.49500(2) 
0.49499(2) 


0.9748(1) 
0.9748(1) 


42.87(20) 


0.0724 


24 


2.1M 


0.38100(18) 


0.38331(14) 
0.38330(14) 
0.38330(16) 


0.9722(7) 
0.9722(5) 


38.0(2.0) 


0.0732 


30 


10. 1M 


0.30645(5) 


0.30765(5) 
0.30764(4) 
0.30757(4) 


0.9701(3) 
0.9700(2) 


32.3(9) 


0.074 


40 


1.4M 


0.21347(12) 


0.21391(10) 

0.21392(9) 

0.21395(8) 


0.9662(9) 
0.9663(5) 


28.2(2.4) 


0.07436 


52 


615fc 


0.15850(16) 


0.15869(13) 
0.15871(12) 
0.15870(9) 


0.9635(15) 
0.9638(9) 


20.5(4.3) 



Table 3: Extracted values of m,m,R, Zn and gn for the Ising model in the symmetric phase. The 
various values for to.r and Zn are given using the :r-space 1-mass fit, naive and momentum space 
fit methods (in this order), 'iter' is the number of sweeps as described in Sect. 4.1 of ref. [11]. 
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K 


L 


iter 


m 


m R 




9r 


5R3 


0.077 


16 


31M 


0.5506(6) 

U.uuUil O I 


0.5589(5) 

U.lJOO t 1 U J 

0.5595(5) 
0.5593(4) 


0.8753(11) 

0.8772(12) 
0.8764(10) 


35.26(12) 


75.7(6) 


0.076 


20 


14.4M 


0.3902(7) 
3890(1 21 


0.3936(6) 
3933(71 
0.3942(5) 
0.3940(4) 


0.8901(20) 
8889(25) 
0.8927(17) 
0.8917(13) 


30.10(13) 


56.2(7) 


0.0754 


32 


6.2M 


0.2606(7) 
0.2594(14) 


0.2619(6) 
0.2616(8) 
0.2622(6) 
0.2621(4) 


0.9007(42) 
0.8988(79) 
0.9028(31) 
0.9020(20) 


25.40(20) 




0.0752 


40 


4.9M 


0.2043(7) 
0.2030(13) 


0.2053(6) 
0.2049(8) 
0.2054(6) 
0.2076(4) 


0.9064(38) 
0.9027(56) 
0.9074(39) 
0.9192(16) 


23.23(24) 


38.8(1.4) 


0.0751 


48 


6.3M 


0.1707(6) 
0.1699(11) 


0.1713(5) 
0.1711(6) 
0.1715(5) 
0.1715(4) 


0.9095(38) 
0.9073(53) 
0.9112(39) 
0.9113(29) 


21.84(22) 


35.7(1.2) 



Table 4: Measured values of m,mB,, Zn, gn and gn 3 for the Ising model in the broken phase. 
Values for m, mR and Z-r, are given using the z-space 1-mass fit, 2-mass fit, naive and momentum 
space fit methods (in this order), 'iter' is the number of sweeps as described in Sect. 4.1 of ref. [11]. 
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